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OUTLINE OF THE TALK 


1- Motivation and Objectives 
2* Formulation 

- Classical and Zero-Total-Pressure-Loss Sets 

- Supersonic Conical Flow Equations 

- Relative Motion in a Rotating Frame of Reference 


3. 

4. 


Highlights of Method of Solution 


Applications : 

- Conical Flow, Sharp-Edged Wings (Classical and ZT PL Sets) 

- Conical Flow, Round-Edged Wings (Classical and ZTPL Sets) 

- Three-Dimensional Flows; Transonic and Low-Speed Flows 

- Uniform Rolling in a Conical Flow 

- Rolling Oscillation in a Locally Conical Flow 


Ex* — p ^ 5 (0 


fs/u*»«r;c mS 


5« Concluding Remarks 



1 Classical 
Vortex 

2 Separation 
Bubble with 
No Shock 

3 No Shock/ 

No Separation 

4 Shock with no 
Separation 

5 Shock-Induced 
Separation 

6 Separation 
Bubble with 
Shock 

7 Vortex with 
Shock 


Figure 1 . Miller and Wood 1 Classification Diagram. 



CLASSICAL EULER EQUATIONS 


• Conservation Form of Euler Equations in a Space-Fixed Frame of Reference 


aq * at * a? . aG n 
U + + Ty + m 0 


(i) 


5 ■ Cp. pu, pv, pw, pe] T 


( 2 ) 


t » [pu, pu 2 ♦ p, puv, puw, puh]* 


( 3 ) 


f • [pv, puv, pv 2 ♦ p, pvw, pvhj 1 


( 4 ) 


5 • [pw, puw, pvw, pw 2 ♦ p, pwh ] 1 


( 5 ) 


e • p/ptr-D ♦ ( u 2 ♦ v 2 ♦ w 2 )/2 


( 6 ) 


h * e + p/p 


( 7 ) 


Conical Variarles 


SUPERSONIC CONICAL FLOW EQUATIONS 


£ • X , n « y/x , C ■ 2 /x 


( 6 ) 


Conical Flow Equations 


si ♦ S * H ♦ 2 e * o 


( 7 ) 


WHERE 


P« F * T) [ 


( 8 ) 


g - G - C E 


( 9 ) 


Zero-Total-Pressure-Loss Euler Equations 


• Replace the energy equation by either 
one of the isentropic gas equations 

p/p Y * const, or + V *(psv) = 0 


Set (1): • Replace the x-momentum equation (second 

elements in the vectors given in Eqs. 
(2) ‘and (3)) by the steady energy 
equation (total constant enthalpy) 

h * const * , r- + \ (u 2 + v 2 + w 2 ) 
TpTTp 2 


Set (2): • Replace the continuity equation (first 

elements in the vectors given in Eqs. 
(2) and (3)) by the steady energy 
equation given in Set (1) 


EXPLANATION OF TOTAL-PRESSURE CHANGE 
FOR CLASSICAL AND ZTPL SETS OF EULER EQUATIONS 

Differential Euler Equations 
Crocco's Theorem 


I ? s * * * ? ♦ ♦ n 

Definition of Entropy Change 


A S * R In - — + C In 

P T P o 


(a) Classicai Set 


Steady flow ^ * o, h - const 


and 


TVS-wxY. AS-Rln 


( 1 ) 


< 2 ) 


For a free-sheet Z is parallel to V! 7 S 3 0 - Pj * Py ♦ Zero-Total-Pressure Loss 
(p) Zero-total-Pressure-Loss Set (Shock-Free and Weak Shocks) 

h * const, 7 S ■ 0 

and must pe parallel to V, Pj * Pj ♦ Zero-Total-Pressure Loss 


Cqhputatqnal Euler Equations 
Crocco's Theorem 


T a s ■ « x v + ^ ♦ vh . I v • ? Viscous-Form of the Equation 


(l) 


Definition of Entropy Change 


a $ » r in * — + cx 


T P n V 


(2) 


(a) Classical Set 

For STEADY FLOW - 0, h » const AND 


TVS»uxV + Numerical Dissipation, A S * R Xn 


Even if w is parallel to V, v S * 0 ♦ Pt * P T ♦ Non-Zero TPL 

* m 1 

(B) Zero-Total- Pressure-Loss Set (Shock-Free and Weak Shocks) 

' h ■ const, A S * 0 

0 » u x V + Numerical Dissipation, P^ ■ P^ * Zero TPL 
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CLASSICAL EULER EQUATIONS FOR THE RELATIVE MOTION IN A ROTATING 
FRAME OF REFERENCE 

• The conservation form of the classical Euler eouations for the absolute motion 

OF THE FLOW IN A SPACE-FIXED FRAME OF REFERENCE 


♦ v • (p v) • o (i) 

? ( P ♦ 7 • (p V V + p T) - 0 (2) 

o t 

H g -l l + V • (p h V) « 0 . (3) 

V 2 

e ■ p/p ( y- 1 ) ♦ y- 

h ■ e ♦ p/p ( 5 ) 


I 



To express these equations in terms of a 
moving frame of reference, we use the following 
relations of the substantial and .local 
derivatives for a scalar "a" and a vector "A": 


Da . D'a 
Ft ' FT" 


aa a d'a 
at at r 



7a 


DA 

Ft 



wxA 


3A a'A 
at “ FT 


V t • 7A + uxA 


(6. a) 


(6.b) 


(7. a) 


(7.b) 


• The transformation velocity V is a function of 

this moving frame of reference translation and 
rotation 


V = V - V = V + wxr (8) 

t r o 


Restricting the motion of the frame of reference to the 
rotational motion, 

OV 

V o = 0 and = q, 

the equations of relative motion in the rotating frame of reference 

+ 7 • (p V ) » 0 
Ot r 

a'(pv) 

§p + 7 * [p V f V p + p I] , _ pfuxf + 2toxV + 


S' (p e ) 

— ap— + 7 • [p h. v J 


where 


v r 1 - - 2 

e r * p{ Y - n * T ~ 1 l uxr i = e - V • (ujxr) 


(9) 


r + o«(wxr)] (10) 

r v r J s - pCV. • (ix?) + (wxr) • (*xF)] (H) 


(12) 


2 

h r ■ + T ■ 7 M ■ h - V • (uxF) 


(13) 




The 

the 


abstract conservative form 
rotating coordinates (x' t y'. 


a 'V, s ' J r »' p r «' 5 r 

-5T^ * -XT' * Sr~ * 


of the relative 
z') is given by 

S 


motion 


in 


terms of 


(14) 


where 


q r - Cp, pu r , pv r . pw r , 

— 2 

£ r * Cpu r . P u r ♦ p, P u r v r . pu r w r . pty.h r J 1 

* Cpv r * pu r v r* pv r + p * pv r V pv r h r jt 
g r * Cpw r . pu r w r . P v r w r . pw p 2 ♦ p, pWr h r ] t 

S * [0. 0. p(J z ♦ anv r + J y] , _ p{ - y + a-v . w 2 

r 

• Since only the rolling motion is solved 

term S has been written for Z = « e , , and u = i e . 

• X ^ 




p(-V r W2 + W f 

the source 


(15) 

(16) 

(17) 

(18) 

• *2 • p - 1 

ojy + ujuy + (JW2 ) J 

(19) 


H16HUGHTS OF METHOD OF SOLUTION 


1. We use the Central-Difference Finite-Volume Scheme with Four-Stage Runge Kutta 

TIME STEPPING AND EXPLICIT SECOND* AND FOURTH-ORDER DISSIPATION TERMS- 

2- For steady flows, local-time stepping is used, and for unsteady flows minimum 

GLOBAL TIME STEPPING IS USED- 

3* A THREE-DIMENSIONAL COMPUTER PROGRAM IS USED TO SOLVE FOR: 

- Conical flows (using 3 conical planes, we enforce the absolute flow vector to be 

EQUAL ON THESE PLANES) 

- Direct solution of the three-dimensional flow problem- 

u. Depending on the problem under consideration, different initial conditions are 
USED- 

5- Depending on the problem ’under consideration, different surface, farfied and 

SYMMETRY COND I T I ONS ARE USED- FOR SUPERSONIC FLOWS, THE OUTER ROW SHOCK IS 
CAPTURED AS PART OF THE SOLUTION- 


NUMERICAL EXAMPLES (1) 


SHARP-EDGED H1NGS (CLASSICAL EULER EOS- ( ZERO-TOTAL-PRESSURE-LOSS SETS) 

ROUND-EDGED MINGS (CLASSICAL EURER EOS- * ZERO-TOTAL-PRESSURE-LOSS SETS) 

- NUMERICAL BOUNDARY CONDITION (COARSE AND FINE 6RIDS) 

- CLOSED FORM BOUNDARY CONDITION (COARSE AND FINE GRIDS) 


THREE-DIMENSIONAL TRANSONIC AND SUBSONIC FLOWS 


original, page ® 

OF POOR QUALITY 


Figure i . Standard Euler Set, Sharp-edged Wing, 64X64 Cell, 

^<»* 2.0, a»10®, (1*70°, £ 2 * 0 . 12 , £^*0.005, 1. Surface pressure, 2 
Crossflow Mach number, 3. Crossflow Velocity. 








Figure 2. Zero-Total -Pressure-Loss Euler, Set (1), Sharp-edged Wing, 64x64 
Cell, M - *2.0, a*10°, p*70°, c 2 *0.12, e 4 *0. 005,1. Surface Pressure 

2. Crossflow Mach number, 3. Crossflow Velocity. 
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Invite Id Shock. 

Induced Separation 




Details ef Laadtng.cdga Flew (r, f . 2.3) 



